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We discuss hadronic light-by-light scattering contribution to the muon anomalous magnetic mo- 
ment aj]^', paying particular attention to the consistent matching between the short- and the long- 
distance behavior of the light-by-light scattering amplitude. We argue that the short-distance 
QCD imposes strong constraints on this amplitude overlooked in previous analyses. We find 
that accounting for these constraints leads to approximately 50% increase in the central value of 
a^', compared to existing estimates. The hadronic light-by-light scattering contribution becomes 
a^i^^ = 136(25) X 10" 
value. 



thereby shifting the Standard Model prediction closer to the experimental 



I. INTRODUCTION 

Recent results I*] from the experiment E821 at BNL 
might indicate a disagreement between the experimen- 
tal value of the muon anomalous magnetic moment — 
{gfi — 2)/2 and the theoretical expectation based on the 
Standard Model (SM). Although no definite conclusion 
is possible at the moment, the experimental value of 
is persistently higher than the SM prediction; the signif- 
icance of the deviation depends on subtle aspects of the 
low-energy hadronic physics. The largest hadronic con- 
tribution to is due to vacuum polarization, see Fig. la. 
It can be found by integrating the (7(e~^e~ —^ hadrons) 
annihilation cross-section with the weight function com- 
puted in perturbation theory. Experimentally, the e^e^ 
annihilation cross section is obtained either from direct 
measurements at low energies or, using the isospin sym- 
metry, from hadronic decays of the r meson. 

The most recent study gives different results for 
the e+e^~ based and the r- based analyses; the primary 
reason is the disagreement between the e'^e~ and the r 
data in the energy range 0.85— 1.0 GeV. It is this experi- 
mental issue that currently limits precision in computing 




FIG. 1: Hadronic contributions represented by quark loops: 
(a) vacuum polarization, (6) light-by-light scattering. 



hadronic vacuum polarization contribution to a^. 

Another source of hadronic contributions to is the 
light-by-light scattering, induced by hadrons, see Fig. lb. 
Compared to the vacuum polarization, this contribution 
is significantly smaller; nevertheless, given experimental 
precision on a^, it is quite important. 

The hadronic light-by-light scattering contribution 
cannot be related to experimental data; for this reason 
the existing estimates of this contribution are model de- 
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pendent. This feature leads to major problems in esti- 
mating both the central value and the theoretical uncer- 
tainty. Given the fact that at low energies the physics of 
light-by-light scattering is non-perturbative, it is nai've to 
expect the fully model-independent solution. The satis- 
factory solution should involve a mixture of both model- 
dependent and first-principles based considerations in 
such a way that the uncertainty caused by the model 
dependence can be minimized and controlled. 

To quantify the quality of the low-energy hadronic 
model, we need a theoretical parameter. Since the per- 
turbation theory is not an option, we must look for the 
parameter other then the QCD coupling constant. The 
two possibilities are the smallness of the chiral symme- 
try breaking and the large number of colors Nc. The 
relevance of these parameters can be seen from the para- 
metrical expression for a 
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ci 



^^QCD 



(1) 



where it is assumed that > m^. Only the power 
dependence on is shown; possible chiral logarithms 
are included into the coefficients ci 2- The first, chirally 
enhanced, term is due to the loops of charged pions in 
the light-by-light scattering. Fig. 2a. The second, Nc- 
enhanced, term is due to exchanges of neutral pion or 
heavier resonances, Fig. 2b. 




FIG. 2: Hadronic contributions to the light-by-light scatter- 
ing: (o) charged pion loop, (fe) exchange of neutral pion and 
other resonances. 

At first sight, it seems natural to expect the chiral pa- 
rameter m^/(47r/7r)^ to be a better expansion parameter 
for a^^^. However, a more careful analysis indicates that 
things can, and perhaps do, work differently. In partic- 
ular, in all hadronic models used to estimate aj]^\ the 
chirally enhanced two-pion contribution is always much 
smaller than the color enhanced contribution. We present 
the "anatomy" of the chirally enhanced 0{N^) contribu- 
tion in the last Section of this paper where we argue that 
this smallness may not be accidental. 

Moreover, a similar example is provided by the 
hadronic vacuum polarization contribution to a^. There, 
the chirally enhanced two-pion contribution gives approx- 
imately 3 X 10~^ which should be compared with the 
iVc-enhanced contribution due to the p-meson that gives 



approximately 50 x 10~^. Although we do not have a 
clear understanding of why the chirally enhanced terms 
are subdominant to such an extent, the above arguments 
suggest that we should accept the dominance of the large- 
Nc expansion over the chiral expansion as the working 
hypothesis. The special feature of the large- A^c limit is 
that scattering amplitudes in any particular channel are 
given by infinite sums of narrow resonances. This helps in 
constructing the model but is clearly insufficient; we need 
further constraints to select among prospective models. 

Such constraints come from the knowledge of short- 
distance behavior of the light-by-light scattering ampli- 
tude, governed by QCD. The asymptotics of this ampli- 
tude at large Euclidean photon momenta is derived from 
the operator product expansion (OPE). The leading term 
in this OPE comes from the quark box diagram enhanced 
by large Nc- This shows a consistency of the OPE con- 
straints with the large-A^c limit. Therefore, we require 
an acceptable large- A"c hadronic model, extrapolated to 
large Euclidean photon momenta, to match the perturba- 
tive light-by-light scattering amplitude. We find that the 
minimal large- A^c model which satisfies this criterion in- 
cludes exchanges of the pseudoscalar 0~ mesons tt^ ,irj,ri' 
and the the pseudovector 1+ resonances oi,/i,/* . It 
is important to emphasize at this point that the model 
with a finite number of resonances is consistent with the 
short-distance constraints for aj]^'; it is known that this 
is not always the case ( see |3| for a recent discussion). 

The short-distance QCD constraints are most restric- 
tive in the pseudoscalar isovector channel. In a special 
kinematic limit, where invariant masses of two virtual 
photons are much larger than the invariant mass of the 
third photon, this channel is completely saturated by the 
neutral pion. The saturation is complete in the sense that 
it works for arbitrary small invariant mass of the third 
virtual photon, in spite of the fact that, in general, the 
OPE applies only when that mass is much larger than 
Aqcd- 

This happens because in the kinematic limit described 
above, the OPE relates hadronic light-by-light scatter- 
ing diagram to the famous "anomalous" triangle dia- 
gram with one axial and two vector currents. Because 
both pcrturbative and non-perturbative corrections to 
the anomalous triangle are absent in the limit of exact 
chiral symmetry, the pion pole contribution is unambigu- 
ous both at small q ~ Aqcd o^nd at large q ^ Aqcd 
momenta. This observation connects the two regions of 
momenta and provides an important constraint thereby. 

In terms of the diagram in Fig. 2a, the constraint 
amounts to the statement that the form factor is present 
in the 7r7*7* vertex if both photons are virtual but it is 
absent if that vertex contains the external magnetic field. 
Although the pseudoscalar channel has been the subject 
of many detailed studies in the past, this constraint has 
been overlooked and, as the result, the 7r°-pole contribu- 
tion to Op was underestimated. This is the main source 
of corrections we find for the pion pole contribution. 

Moreover, additional constraints on subleading terms 
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in 7r7*7* form factor, derived long ago in Ref. 4] , were not 
utilized previously. Accounting for these constraints, also 
leads to the increase in the result. As a consequence, the 
central value of the pion pole contribution to increases 
by approximately 20 x 10^^^. Similar increases occur 
for other pseudoscalar [rj,!]') and pseudovector channels 
(ai,/i,/r). 

Unfortunately, the constraints on all, but tt^, ex- 
changes are not very restrictive; because of that we can- 
not claim significant reduction in the theoretical uncer- 
tainty of hadronic light-by-light scattering contribution 
to a^. Nevertheless, imposing all the constraints from 
the short-distance QCD, we arrive at aj]^' = 136(25) x~^^ 
which is approximately 50 per cent larger than the exist- 
ing estimates '5','^[1|||. 

The rest of the paper is organized as follows. In the 
next Section we discuss the constraints coming from the 
short-distance QCD and a minimal model for hadronic 
contributions to a^^^. We consider the pseudoscalar and 
the pseudovector exchanges in Sections III and IV, re- 
spectively. In Section V we briefly discuss the C(iV°) 
pion box contribution to aj}''. We present our conclusions 
in Section VI. Additional formulas are given in Appen- 



dices. 



II. SHORT-DISTANCE QCD CONSTRAINTS 
AND HADRONIC MODEL 

In this Section we describe the constraints coming 
from the short-distance QCD and formulate the hadronic 
model that satisfies these constraints. 



A. Kinematics 

We begin with the kinematics. The light-by- light scat- 
tering amplitude involves four photons with momenta qi 
and the polarization vectors e^. We take the photon mo- 
menta to be incoming, '^qi — 0. The first three photons 
are virtual, while the fourth one represents the external 
magnetic field and can be regarded as a real photon with 
the vanishingly small momentum q^. The amplitude M 
is defined as 



= -e'l d'^x d^y e-«^--«=« e^,^ 6^ ef (G|T {j,, (x) j,, (y) j,, (0)} I7) , (2) 

I 



where is the hadronic electromagnetic current, = 
qQ"/^q, written in terms of the three quark flavors q — 
{u, d, s} with Q being the 3x3 diagonal matrix of quark 
electric charges. In addition, f^^ = qje^ — q^ej denotes 
the field strength tensor of the soft photon; the light-by- 
light scattering amplitude is proportional to this tensor 
due to gauge invariance. Since A4 is linear in the small 
momentum 94, for the purpose of computing the light-by- 
light scattering contribution to a^x, we can set 54 = in 
the tensor amplitude -^^i/iaMaT* calculate it assuming 
that (71 -f 92 + 93 = for the virtual photons. Because the 
momenta 91,92,93 form a triangle, there are just three 
independent Lorentz invariant variables; we choose them 
to be the virtualities of the photons 94_3. 

In general, the light-by-light scattering amplitude is 
a complicated function of photon's virtualities. How- 
ever, there are only two distinct kinematic regimes in 
the light-by- light scattering amplitudes: the Euclidean 
momenta of the three photons are comparable in mag- 
nitude 9i ~ 92 ^ 93, or one of the momenta is much 
smaller than the other two. The second limit can be 
analyzed in a very simple fashion using the OPE of the 
light-by-light scattering. Also, this limit is of importance 
because it helps us to identify the pole-like structures in 



the OPE amplitudes and in this way connect the OPE 
to phenomenological models. 

B. OPE and triangle amplitude 

Since the light-by-light scattering amplitude is sym- 
metric with respect to photon permutations, we can 
study the second limit assuming that 9i « 92 ^ 93 • 
In this kinematic regime, we begin with the well-known 
OPE (see e.g. 0) for the product of two electromagnetic 
currents that carry the largest momenta 91, 92, 

I J d^x dS e-^^--n^v T {j^^ (y)} = 

y'd4ze-(^^+^^)^|^6^,^,,pg^j|(z)-|-... . (3) 

Here, — 9 (3^7'' 75 9 is the axial current, where different 
flavors enter with weights proportional to squares of their 
electric charges and 9 = (91 — 92)/2 w 91 w —92 . We 
retain only the leading (in the limit of large Euclidean q) 
term in the OPE associated with the axial current ; the 
ellipsis in Eq.(|3Jl stands for subleading terms suppressed 
by powers of Aqcd/9- The momentum 91+92 = —93 
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flowing through is assumed to be much smaller than 
q. We note in passing that Eq.© has been apphed earlier 
in various situations; for example, the matrix element of 
Eq. (|3J) between the pion and the vacuum states gives the 
asymptotic behavior of the 7r'^7*7* amplitude at large 
photon virtualities Q. 

For the purpose of further discussion it is convenient 
to present the current j^p as a linear combination of the 

isovector, j^^^ = qXa^fp-fsq, hypercharge, jj^^^ = qXajpjsq, 

and the SU(3) singlet, — qjpj^q, currents, 

. _ Tr [XgQ^] .(,) 

^^p- 2^ Tr[A2] -^sp 



a=3,8,0 



(4) 



where Ao is the unity matrix. 

Once the dependence on the largest momenta qi 2 is 
factored out, the next step is to find the dependence of 
the light-by-light scattering amplitude on the momentum 
qs. This dependence is given by the amplitudes T-^p^ 
that involve axial currents j^'p' and two electromagnetic 
currents, one with momentum (73 and the other one (the 
external magnetic field) with the vanishing momentum 



T{jii{z)jpM}\l)- (5) 



The triangle amplitudes for such kinematics were con- 
sidered recently in 11] . It is shown in that reference that 

T^p^ can be written through two independent amplitudes. 



rp{a) 

M3P 



-^2 Y"l \%) <l3pq3Ut^a + 

rKqD^-qlfpsp+qspsisUp-qspqsU,.,)]- (6) 



The first (second) amplitude is related to the longitudinal 
(transversal) part of the axial current, respectively. In 
terms of hadrons, the invariant function wl(t) describes 
the exchanges of the pseudoscalar (pseudovector) mesons. 

In perturbation theory wl.t are defined by the famous 
triangle diagram. For massless quarks, we obtain: 



wl'{q ) 



2w'-^\q^) 



q- 



(7) 



An appearance of the longitudinal part is the signa- 
ture of the axial Adler-Bell-Jackiw (ABJ) anomaly [lOl |. 
Although the perturbation theory is only reliable for 
q^ ^ ^QCD' where it coincides with the leading term 
of the OPE for the time-ordered product of the axial and 
electromagnetic currents, the expressions for longitudinal 
functions w^^'^^ given in Eq.Q are exact QCD results in 
the chiral limit ruq = for nonsinglet axial currents. The 
fact that there are no perturbative 12] and nonpertur- 
bative [T^ corrections to the axial anomaly implies that 

in Eq.jTJl is correct all the 



the pole behavior of wjf'*' 

way down to small q^, where the poles are associated 



with Goldstone pseudoscalar mesons, vr" in w^^^ and rj in 

(8) 
Wl ■ 

Equations 10 and iQ allow us to derive the coupling 
of the tt" meson to photons. To this end, consider the 
isovector part of the triangle amplitude T^p . The residue 
at q^ = 0, corresponding to the tt" pole, is the product 
of two matrix elements, 

(Oljf^lTr") - 2^F^qp , (7r"ljp3l7) = -4e5.^^g^/.;.3 -(8) 

Comparing with Eqs.lQ, we derive the well-known 
result |1Q] for 7r77 coupling: 



^cTr [AaQ^ 
167r2 



(9) 



In a similar way, the g-ri-y-y coupling in the chiral limit can 
be derived, if needed. 

The absence of perturbative and non-perturbative cor- 
rections and therefore the possibility to use the OPE 
expressions for vanishing values of q^ is unique for the 
longitudinal part of nonsinglet axial currents.^ For the 
transversal functions wt as well as for the singlet longi- 



tudinal function 



(0) 



there are higher order terms in the 



OPE that, upon summation, generate mass terms that 



shift the pole position 1 /q- 1/(5 



We use this 



modification of the pole-like terms for each channel in 
what follows. The lightest pseudovector mesons are the 
ai(1260), /i(1290) and A* (1420) mesons. For the singlet 
axial current the pole in wf ^ is shifted to m^, 

Consider a triangle amplitude for any isospin channel 
in the limit of large q^, where the OPE and the per- 
turbation theory are applicable and Eq.jTI) is valid. An 
important consequence of this equation is that triangle 
amplitudes are not suppressed for such values of q^. In 
terms of hadrons it means that no form factor is present 
in the ^17*7 interaction vertex where the real photon is 
soft (external magnetic field). This is in clear contradic- 
tion with the common practice 0, 0) when, for tt^ ex- 
change, the form factor FT^~ft^(q'^ ,0) is introduced. Such 
transition form factor has to be present when one of the 
photons is virtual, the other photon is on the mass shell 
and the pion is on the mass shell as well. However, this is 
not the kinematics that corresponds to the triangle and 
the light-by-light scattering amplitudes, relevant for 
computation, where the pion virtuality is the same as 
the virtuality of one of the photons. The absence of the 
suppression is also consistent with the dispersion repre- 
sentation of the amplitude since the imaginary part is 
nonvanishing only at q^ = in the chiral limit. 



^ More precisely, pert urbative corrections to Wj, 
as shown in Ref. ll4l . 



(3,8) 



are also absent 
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Combining Eas. (|3l6(l . we write the the hght-by-hght ampHtude -4pi;^2A'37i5 fo'" Qi ~ oi ^ 1z following form: 

/^P3 + +Q3^s<i3fS-<ik3Us) } + •••, (10) 



a=3,8,0 



r 



where no hierarchy between and Aq(-,j-, is assumed. 
The weig hts PF^") are defined as 



(Tr [A,Q2 
Tr[A2]Tr[Q4 



(11) 



4' 12' 3 



In the limit g| ^ ^qcD' Ea. HlOl) can be simplified us- 
ing the asymptotic expressions Eq.(jZIl for the invariant 
Convoluting the tensor amplitude with 



functions 



(a) 

"l.t- 



the photon polarization vectors and analytically contin- 
uing to Euclidean space, we arrive at: 



4 . . 4 / 

-2-^ { /2 /l } { //3 } - -2"^ ({ (72 /2 /l //a 93 } 
q^q q^H ^ 



+{ 91/1/2 //3 93}- 



9? 



{/2/1K//3} 



Here, f^'^ = gfej' — q^e^ are the field strength tensors, 
the braces denote either traces of products of the matrices 
/f or their convolutions with vectors qt. 

In Ea. (|12|l and in the remainder of the paper, we use 
Euclidean notations instead of Minkowski ones used be- 



fore. The continuation to Euclidean space mostly con- 
cerns the change in sign for all qf and the overall change 
in sign for the amplitude A, since it involves the product 
of two Levi-Cevita tensors. The result can be verified 
by comparison with the direct computation of the quark 
box diagram, for arbitrary qi_^, presented in Appendix 
I. There we show that the amplitude can be described 
in terms of nineteen independent tensor structures and 
five independent form-factors. In what follows, we mostly 
deal with the approximate form of the amplitude Ea. (|12|l . 
but we make occasional references to general expression 
in Appendix I. 



C. The model 

Two different terms in Ea. (|10|l can be identified with 
exchanges of the pseudoscalar (pseudovector) mesons 

(12) for the functions w'f'lplq'^) . Extrapolating Ea. (|12|l from 

Qi 2^ Aqqj) to arbitrary qf 2, we arrive at the following 



model: 



where 



A = Aps + ApY + permutations. 



(13) 



J 



a=3,8,0 



^PV - W^^"V^"^('Z?,g|)4"^('Z3') ({92/2/l/73<Z3} + {gi/l/2/73<Z3} + 

a=3,8,0 ^ 



{f2fl}{fh} 



(14) 
(15) 



r 



The form factors 02^"r(9i, f?!) account for the dependence 
of the amplitude on g^2- Pictorially (see Fig. 2b), these 
form factors can be associated with the interaction vertex 
for the two virtual photons on the left hand side, whereas 
the meson propagator and the interaction vertex on the 
right hand side form the triangle amplitude described by 
the functions w^^lp^q"^). In the next Sections we intro- 
duce models for these functions consistent with the short 
distance behavior of the light- by- light scattering ampli- 
tude. 

Note that our model does not include explicit ex- 
changes of vector or scalar mesons. This is a consequence 
of the fact that, to leading order, the OPE of the two 



vector currents produces the axial vector current only. 
However, the vector mesons are present in our model im- 
plicitly, through the momentum dependence of the form 



factors 



as well as the transversal functions 



III. CONSTRAINTS ON THE PSEUDOSCALAR 
EXCHANGE 



The tt" exchange provides the largest fraction of the 
hadronic light-by-light scattering contribution to a^. It 
is therefore appropriate to scrutinize this contribution 
as much as possible and ensure that it satisfies all the 
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possible constraints that follow from first principles. 

As we discussed earlier, the longitudinal part of the tri- 
angle amplitude is fixed by the ABJ anomaly. Account- 
ing for explicit violation of the chiral symmetry given by 
the small mass of the pion, we derive 



The ABJ anomaly also fixes 0^^'' (0, 0) 



•/-f (0,0) = 



47r2F2 ' 



(16) 



(17) 



which is the correct asymptotics indeed [J. This means 
that the neutral pion exchange in Ea. H14ll saturates the 
corresponding short- distance QCD constraint. 



This comparison also proves our previous claim that 
the form factor _F!n-^*^(g|, 0) cannot be present in the 
amplitude Ea. ((TH|l : if that form factor is introduced, the 
asymptotics of the light-by-light scattering amplitude be- 
comes l/gl, as opposed to l/gf behavior that follows 
from perturbative QCD. This proof is, of course, equiva- 
lent to our discussion of the triangle amplitude in Section 
11. 



so that the model for the pion exchange in the light-by- 
light scattering amplitude takes the form. 



^71-0 = 



27r2F2 g2 _^ ^2 
' permutations . 



{/2.fl}{/73} 



The 7r7*7* form factor i^7r7*7' (9i , ^l) is defined as 



2 



'/'^(O-O) 



(18) 



(19) 



The comparison with the OPE constraint given by the 
relevant term in Ea. lfTUIl leads to 



lirn F^^.^.{q^ ,q^) 



8^2 



(20) 



The absence of the second form factor in the amplitude 
Ea. (|18(l distinguishes our approach from all other calcu- 
lations of the pion pole contribution to that exist in 
the literature. As we show below, it has a non-negligible 
impact on the final numerical result for the pseudoscalar 
contribution to a^. Here we note, that the result for the 
pion pole contribution is expected to increase, because 
the absence of the second form factor leads to slower 
convergence of the integrals over loop momenta, making 
the result larger. As we show below, this is indeed what 
happens. 



Further constraints on the model follow from restric- 
tions on the pion transition form factor Fj^jty* that were 
recently reviewed in 8]. For numerical estimates we use 
their LMD+V form factor 



47r2F2 qfqliqf 



4) 



h2qlql + hM + «i) + (A^cM4M|/47r2F2) 



[ql + Ml){ql + Mj){ql + Ml){ql + Mj) 



(21) 



r 



where Mi = 769 MeV, A'h = 1465 MeV, /is = 
6.93 GeV''. 

The parameter /12 was not determined in Ref. and 
we can fix it if we notice that it contributes to the l/q^ 
correction to the leading asymptotics of the pion form 
factor, Ea. H2(J|l . Such correction comes from the twist 4 
operators in the OPE expansion of the two electromag- 
netic currents Eq.©- It was analyzed long ago in Ref.i^] 
using the OPE and the QCD sum rules approaches. The 
result of such an analysis implies that the coefficient of 
the 0{q~'^) term in the asymptotics of the pion form fac- 
tor is numerically small; in terms of the parametrization 
Ea. (|21|l . this means that /12 ~ —10 GeV^ has to be cho- 
sen. We use this value for numerical estimates in what 
follows. 

Equations H18|) and (|21|l completely specify the model 
for the pion pole contribution that we use for numerical 
calculations below. Before going into that, we discuss the 



sensitivity of the final result to possible modifications of 
the model. 

We denote the structure that multiplies {/2/i}{//3} 
in Eq.lHl) as W^^^G^°'^{ql,ql,ql). Comparing the 7r°- 
pole exchange amplitude, Ea. H18|l . to the full light-by- 
light scattering amplitude (see Appendix I), we find that, 
for asymptotically large virtualities of the photons, the 
matching requires 



GT%lql<ll)=G2{qlqlql). 



(22) 



Consider Eq.(|221l in the hmit Aq(-.j3 < (Zi < (?| ~ <?|. 
It is easy to sec that the left hand side in Eq. (|22ll develops 
the I/qi behavior; from expression for G2 in Appendix 
I it follows that G2{q2,q2TQi) ~ l/^i such kinematic 
regime. Hence, there is a mismatch between our model 
and the OPE prediction. 

The second option is to consider Ea. (|22|l in the situa- 
tion when all the momenta are asymptotically large and 
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equal in magnitude (Zi = = 93 = 9^ • this regime 

8 



3g4 



whereas 



(23) 



(24) 



Again, the model fails to describe the OPE constraint 
perfectly. 

Of course, the above failures do not necessary invali- 
date the model; after all we are interested in the light- 
by-light scattering contribution to and various regions 
of loop momenta contribute differently to the integral. 
For this reason, we have to investigate if the above mis- 
matches influence the numerical estimate for the pion 
pole contribution to aj]^^ To this end, we notice that 
G'™°'^(g|, gfj ^i) can be modified by adding to it 



factor.^ The 'r]{r]')j*j* interaction vertex is normalized 
in such a way that the decay widths of these mesons into 
two photons is correctly reproduced; this allows to ac- 
count for the rj — 1]' mixing in a simple way. 

How good these "experimental" couplings compare to 
the theoretical expectations based on our model? Be- 
cause of the rj — rj' mixing, we expect that the sum of 
viv')!! couplings squared is predicted by the model more 
accurately than each of the couplings separately. We find 



^7777 



= 3, 



(26) 



SCr^iqlqlql) 



^qhl 



whereas using experimental values for the r]{r]')"f"f cou- 
plings we arrive at r = 2.5(1). Although we use this 20% 
discrepancy as an error estimate on the rj + rj' contribu- 
tion, we note that it rather implies an increase in the 
result since the agreement between "experimental" and 
theoretical asymptotics can be improved by adding more 
pseudoscalar mesons to the model. 

Compared to the results quoted in 8] , removal of the 
second form factor increases the ry and rj' contributions 



{qf + Mt){qi + Mi){qi + M()iqi + Mi) ' 



(25) from approximately 



a 



n 



13 X 10-" to a2 



a]] = 18 X lO"-"^^. The sum of the contributions from all 



without running into a contradiction with the required pseudoscalar mesons {tt° ,7],!]') leads to the estimate: 
pole behavior with respect to q^. After adding SG™°'^, it 



is easy to see that, by tuning ^, one can either ensure that 
the pole in qf is absent or that the asymptotic behavior 
of 0^°"^ + (5G^°^ becomes consistent with Eq.^. The 
two constraints are satisfied for ^ = — 2 and ^ = —4/3, 
respectively. 

We can investigate the importance of these constraints 
by computing the contribution of dG'^°'^ to aj]^' for ^ = 1. 
Upon doing so, we find that it changes by « 0.5 x 
10""'^^. Hence, regardless of the value of ^, (5G™°'* can be 
neglected at the current level of precision. We therefore 
use Eas. H18l21|) as our model for the pion form factor in 
what follows. 

The result for a'^ with the LMD+V form factor for 

/i2 = -10 GeV^ quoted in 8] is af = 63 x 10"". Using 
the formulas in 8] it is easy to repeat their calculation 
removing the pion transition form factor that involves 
the soft photon. In that case, the result becomes 76.5 x 
10"^^, a shift in the positive direction. In addition, as we 
mentioned earlier, we consider the value /12 = —10 GeV^ 
to be preferable because of the OPE constraints on the 
pion transition form factor. Note, however, that h2 = 
was used in Q to derive the central value 58 x 10~^^; 
compared to that number, our central value is larger by 
approximately 20 x 10~^^. 

A similar analysis for the isosinglet channels leads to 
the conclusion that these channels are saturated by 77 
and 77' mesons; matching to pQCD result suggests that 
no transition form factor is present for the soft photon 
interaction vertex in those cases as well. Since these con- 
tributions are smaller than that of ttq, we do not use 
sophisticated models for i] and rj' transition form factors 
and estimate them using the simplest possible VMD form 



flPS = 114(10) X 10" 



(27) 

The central value in Eq. (|27|l is almost 40 percent larger 
than most of the existing results for a' 

IE 110 The 

major effect comes from removing the form factor for 
the interaction of the soft photon (magnetic field) with 
the pseudoscalar meson; the necessity to do that unam- 
biguously follows from matching the pseudoscalar pole 
amplitude to the pQCD expression for the light- by- light 
scattering. 

On the contrary, the error estimate in Ea. H27() is sub- 
jective; it is based on the variation of the result when 
input parameters of the model are varied. It is impossi- 
ble to defend the exact number for the error estimate in 
Ea. (|27() : however, we believe that it adequately describes 
our current knowledge of the pseudoscalar contribution. 



IV. PSEUDOVECTOR EXCHANGE 

In this Section we discuss the pseudovector exchange 
amplitude Apy, Eq.lO. From Eqs. Q, lO, we 

find the asymptotics of 0^'' and w^'' , 



lim 



-4 



lim 



(a) 



iq' 



(28) 



^ The VMD form factor obviously violates the scaling of the form 
factor when both photon virtualities become large. We have 
checked that using the form factor consistent with the asymptotic 
scaling at large values of gi,2 has no bearing on the final 

result for the r] and rj' contributions. 
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As we mentioned earlier, the lightest pseudovector reso- 
nances are the Oi meson with the mass M^^ = 1260 MeV, 
the /i meson, with the mass Mf-^ = 1285 MeV and the /j* 
meson with the mass Mf* = 1420 MeV. The contribu- 
tion of these mesons to al!^' is cut off at the scales defined 
by their masses. This suggests that the pole-like singu- 
larities in Ea. H12|) are shifted from zero to the masses of 
the corresponding pseudovector and vector mesons. We 
also remind the reader that (jj^^qf, describes the form 
factor for the 7*7* transition. Shifting all the 

poles by the same amount, i.e., neglecting mass differ- 
ences, we get the simplest possible model consistent with 
perturbative QCD constraints Ea. H28|l . 



(a) 



(9) 



(29) 



This implies, in particular, that we do not distinguish 
between different isospin channels. 

Although this model is not very realistic, we can use 
it to derive a simple analytic result which will help us 
to exhibit the dependence on the mass scale M. Assum- 
ing that M 3> mfj_, we compute the contribution of the 
pseudovector meson to and obtain: 



71 

192^ 



81. 

le*^' 



144 



Also, we note that the result Ea. H30|l exhibits strong 
sensitivity to the mass of the pseudovector meson and the 
mass parameter in the form factor. If we associate the 
mass scale M in Ea. (|30|l with the mass of the p-meson, 
the result increases roughly by a factor 4 and becomes 
flP^ ~ 28 X 10~^. Because of the strong sensitivity to 
the mass parameter, we have to introduce a more so- 
phisticated model accounting for the mass differences in 
different isospin channels. 



Let us start with the isovector function w, 



(3) 



This 



function describes the triangle amplitude that involves 
the isovector axial current, the virtual photon and the 
soft photon. We expect therefore that w^\q'^) should 
contain two poles with respect to gf: the first one, that 
corresponds to the ai(1260) pseudovector meson and 
the second one, that corresponds to the vector mesons 
p, w, thereby reflecting the properties of the virtual pho- 
ton. Such a model was constructed in Ref.[Tl| where 
it was required that, for large values of q^, the equality 
Whiq^) — 2wt(9^), remains valid through 0{q~'^) terms. 
Such a requirement leads to 



(31) 



where we do not distinguish between the masses of p and 
Lj mesons. Correspondingly, the form factor (/'^^■" (gj, gf) 
becomes 



TfL 

1010 ^ X 10' 
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(30) 



ml) 



(32) 



where S2 = 0.26043. Using M = 1300 GeV as an exam- 
ple, we obtain aF^ = 7 x 10"-^^. 

There are two comments we would like to make about 
this result. First, we compare it to the existing estimates 
of the pseudovector meson contribution |^ U- In those 
references, the results 2.5 x 10^^^ and 1.7 x 10~^^ have 
been obtained. We have checked that the difference be- 
tween our result Ea. (|30|l and the results of 5, 7] can be 
explained by the absence of the form factor for the 7*7/1 
interaction vertex in our model; when such a form factor 
is introduced, our result decreases to 2.6 x 10~^^, in good 
agreement with the estimates in |^ |^ . 



For the isoscalar pseudovector mesons /i(1285) and 
/i(1420) we assume the "ideal" mixing similar to uj and 
(/); this assumption is consistent with experimental data 
for decays of these resonances. Then, instead of the hy- 
percharge and the SU(3) singlet weig hts VK(*) and W^°\ 
we use 



36' 18' 



(33) 



and the following expressions for the corresponding func- 
tions wt and 4>t ■ 



w. 



{u+d) 



v^\q') = — 



"^/i " (™y5) ml - (m2/5) 



ml -I- m. 



q^ +ml 



i(u+d) / -i 2\ 



iql+ml){q. 
4 



-ml){q^+ml) 



ml)' 



(34) 



Note, that these refinements of the simple expression for 



the function wt in Ea. l|^ make the effective mass of 
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the pseudovector meson lower. This leads to the increase 
in as compared to Ea. (|30|l . We obtain the following 
estimate: 

flPV ^ 7 ^ Q _^ 0.8) X 10"" 22 X 10"", (35) 

where the three terms displayed separately are due to the 
isovector, u + d and s exchanges respectively. 

To check the stability of the model, we consider an 
opposite case for the mixing, assuming that /i is a pure 
octet and and /j* is an SU(3) singlet meson. The estimate 
for then becomes 

^ 7 ^ i Q _^ 9 7) X 10-" = 17 X 10"". (36) 

We see that the SU(3)-singlet contribution is significant, 
in spite of the fact that the corresponding masses are 
the largest. The reason for such a behavior is a stronger 
coupling of the SU(3)-singlet meson to two photons. We 
see also that in spite of a very strong redistribution be- 
tween the different SU(3) channels, the final result for 
the pseudovector contribution is relatively stable against 
such variations of the model. 

We use the result for the pseudovector contribution in 
Ea. (|35|l in our final estimate of aj]^' assigning ±5 x 10^^"'^ 
as an error estimate. 



V. THE ANATOMY OF THE PION BOX 
CONTRIBUTION 

In this Section we make a few comments concerning 
another contribution to aj}'' frequently considered in the 
literature, the so-called pion box contribution. This con- 
tribution is peculiar because, being independent of the 
number of colors Nc, it is enhanced by the other po- 
tentially large parameter, the small value of the pion 
mass relative to the scale of chiral symmetry breaking 
~ 1 GeV. 

The results for the pion box contribution to aj]^' were 
obtained in H 0; they are aP'°" = -4.5(8.5) x 10"" 
in and oP'™ = -19(5) x 10"" in 0. The differ- 
ence between the two results is attributed to a different 
treatment of subleading terms in the chiral expansion; 
while the vector meson dominance (VMD) model is used 
in to couple photons to pions, the so-called hidden 
local symmetry (HLS) model is used in Although 
the smallness of aP'°" shows that the chiral enhancement 
is not efhcient for ajjj'', the strong sensitivity of the final 
result to the particular method of including heavier reso- 
nances suggests that the chiral expansion per se may not 



The claim in m and H that the VMD model violates the Ward 
identities for the 7*7*7777 amplitude is not correct, if the VMD 
is implemented in the standard way, by introducing the factor 
{M^Qixi/ + g^(j'i^)/(M^ + q^) for each photon in any interaction 
vertex. The Ward identities, discussed in Q, are then automat- 
ically satisfied. 



be a reliable tool for this problem. If this is true, the 
natural question is to what extent the above estimates of 
the pion box can be trusted. With this question in mind, 
we investigated an "anatomy" of this contribution based 
on the analytic calculation of aP'°" in the framework of 
the HLS model. 

The logic which is behind the use of the chiral expan- 
sion to estimate subleading 0{N^) contributions to a}}'' 
is as follows. If the pion box contribution to is deter- 
mined by small values of virtual momenta, comparable 
to the masses of muon and pion, we can compute it by 
using chiral perturbation theory. The leading term in the 
chiral expansion delivers a parametrically enhanced con- 
tribution (a/7r)^(m^/m7r)^ to a™, which can be derived 
from the scalar QED Lagrangian for the pions: 

Df,^:D^'^r* -ml\7r\^. (37) 

Here Z?^ = 5^ — ieAf^ is the covariant derivative and tt 
is the pion field. The Lagrangian Ea. H37|) is the lead- 
ing term in the effective chiral Lagrangian and hence 
the terms neglected in Ea. l|T7|) . are suppressed by the 
square of the ratio of the pion mass to the scale of the 
chiral symmetry breaking. Numerically, these correc- 
tions are expected to be small since m^/M^ ^ 0.04 and 
mj/(47r/^)^ ~ 0.025; therefore, they should not change 
the scalar QED prediction by more than a few per cent. 

It is then puzzling that the results available in the lit- 
erature exhibit drastically different behavior. Existing 
calculations show that the scalar QED contribution is re- 
duced by a factor from three to ten , when sublead- 
ing terms in the chiral expansion are included. Hence, 
the results for the pion box contributions existing in the 
literature tell us that the chiral expansion for this con- 
tribution does not work. In order to identify the reason 
for that, we computed several terms of the expansion in 
rriTr/Mp in the framework of the HLS model. Comparing 
the magnitude of the subsequent terms in the expansion, 
we can determine the rate of convergence of the chiral 
expansion and estimate the typical virtual momentum in 
the pion box diagram. 

As we demonstrate below, the typical virtualities in the 
pion box diagram are approximately 4m7r which leads to 
a slow convergence of the chiral expansion and explains, 
to a certain extent, a very strong cancellation between 
the leading order scalar QED result and the first /Mp 
correction. The remaining terms in the chiral expansion 
are smaller (although not negligible). 

Large value of typical virtualities brings in another 
problem with the scalar QED model Eq. H37|) and its mod- 
ifications based on the VMD. Since fairly large virtuali- 
ties are involved, one might wonder about the quality of 
the model for asymptotically large values of q. To see 
that the model fails relatively early, we can consider the 
deep inelastic scattering of a virtual photon with large 
value of q^, on a pion. The Lagrangian (|37|l then implies 
the dominance of the longitudinal structure function, 
while QCD predicts the opposite. Modifying the scalar 
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QED Lagrangian Ea. (|37|l to accommodate the VMD ei- 
ther directly or through the HLS model, does not fix this 
problem, since only an overall factor {Mp/{Mp + q^)Y is 
introduced in the imaginary part of the forward scatter- 
ing amplitude. This mismatch implies that the models 
that we use to compute the pion box contribution be- 
come unreliable very rapidly once the energy scale of the 
order of the p-meson mass is passed. Since 4m7r is only 
marginally smaller than Mp, it is hard to tell how big a 
mistake we make by ignoring the fact that our hadronic 
model has incorrect asymptotic behavior. 

The above considerations suggest that while it is most 
likely that the pion box contribution to is relatively 
small, as follows from a strong cancellation of the two 
first terms in the chiral expansion, the precise value of 
this contribution is impossible to obtain, using simple 
VMD and the like models. 

We now perform an analytic calculation of the pion box 
contribution to and demonstrate that the typical loop 
momenta in the pion box amplitude are relatively large. 
For the analytic calculation, we use the HLS model 5] to 
describe low-energy hadron- photon interactions. From a 
computational point of view, we have to deal with three- 
loop diagrams that involve three distinct scales: the mass 
of the muon m^^ the mass of the pion rriT^ and the mass 
of the p meson Mp. Because the masses of the muon and 
the pion are close, one can treat them as almost equal 
and construct an expansion in their mass difference; this 
reduces the problem to two-scale diagrams. As the next 
step, one constructs the expansion in m^/Mp, using the 
theory of asymptotic expansions for Feynman diagrams 
(see jial for a review). 

As it turns out, there are twelve different momenta re- 
gions to be considered; the two limiting cases are a) all 
the loop momenta are much smaller than the mass of the 
p-meson and b) all the loop momenta are comparable to 
the mass of the p meson. In case a) one has to compute 
the three-loop "on the mass-shell" diagrams; in case b) 
the masses of both muon and pion can be neglected and 
one has to compute the three loop vacuum bubble di- 
agrams. Intermediate cases in which some of the loop 
momenta are small and the other are large factorize into 
the product of one- and two-loop diagrams. The tech- 
niques needed for such a computation are described in 
Refs.[ll[l3. 

We now present the result of the calculation. To do 
this in a compact form, we introduce the notation 5 ~ 
(m^ — mTj)/mTj and L = \n{M p / m.^) . We then write: 

i=0 ^ P 

The functions fi{5, L) for i = 0, 1, 2 are given in Appendix 
II. We have computed fi{6,L) for i from j = to « = 4 
analytically and we use those functions below for numer- 
ical estimates. In addition, we use = 136.98 MeV, 
= 105.66 MeV and Mp = 769 MeV. With these input 



values, Ea. (|38|l evaluates to: 

= -0.0058 ^ = (-46.37 + 35.46 + 10.98 
-4.70 - 0.3 + ...) X 10"" = -4.9(3) x 10"", (39) 

where the subsequent terms in Eq. (|39|l correspond to the 
subsequent terms in Ea. H38|) . 

The feature of the result Ea. (|39|) which has to be 
noticed is the strong cancellation between the leading 
and the subleading terms in the chiral expansion; the 
other terms, being non-negligible numerically, are cer- 
tainly smaller. It is this cancellation that ensures the 
smallness of the final value for the pion pole contribution 
to . We can use Eq. I|39() to determine typical momenta 
virtualities in the pion box contribution. 

For simplicity, we study this question assuming m-rr — 
m^, which implies 5 = in the formulas for /i(<5, V) pre- 
sented in Appendix II. In this limit Ea. H39|l becomes: 

a1{m^ = nip) « (-69 + 54 + 18 - 8 - 1 + ..) x 10"". 

(40) 

We also assume that the contribution to can be de- 
scribed by the chiral expansion, with the effective scale 
/i. This scale characterizes the typical virtual momen- 
tum in the pion box diagram. Motivated by the chiral 
perturbation theory, we make an Anzats: 

a,(™. = m,) « (-j — [c, + c— + c— + ... j . 

(41) 

We further assume that all the coefficients in the above 
series are numbers of order one. Setting ci = 1 in the 
above equation, we can determine the value of fi by 
comparing it with the first term in Ea. l40|) . We obtain 
/i = 4.25™^. Then, Ea. (|41|l becomes: 

al{m^ = mp) « (-69 + 41c2 + 24c3 + 14c4 + ..)x 10"", 

(42) 

which implies that with C2 ~ 1.3, C3 ~ 0.75 and C4 « 
—0.6, we can easily fit Ea. (|^ . 

The above calculation suggests a simple way to under- 
stand the magnitude of the the chirally suppressed terms 
in Eq.||211). Since p. « 4m^ « 550 MeV < Mp, the chiral 
expansion converges, but rather slowly. Therefore, the 
estimates based on the chiral expansion do make sense 
in principle. A closer look at fi {6, L) reveals that these 
functions contain In (Mp/m7r) -enhanced terms. However, 
in view of the above argument, the appropriate way to 
write the large logarithms is ln{Mp/p); doing so, we ob- 
serve that "large" logarithms become rather moderate 
numerically and every function fi (5, L) is dominated by 
constant terms. 

We therefore see that the typical virtual momenta in 
the pion box contribution are larger than the mass of 
the pion by, approximately, a factor of 4. While the chi- 
ral expansion is still a valid tool for such virtualities, its 
predictive power becomes small. This can be seen from 
Ea. (|39|) . which implies that the final result for the pion 
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box contribution to ajy is very sensitive to higher order 
power corrections. It is clear that since none of the mod- 
els, be it the HLS model or the VDM model, can claim 
full control over higher order power corrections in the 
chiral expansion, the exact result for the pion box con- 
tribution is not very meaningful. However, the fact that 
the chiral expansion is still applicable suggests that the 
strong cancellation between the leading order term and 
the first subleading term in the chiral expansion may be 
a generic feature of QCD. 

Therefore, we find it reasonable to believe that the pion 
box contribution to a!!'' is much smaller than the estimate 
based on the chirally enhanced scalar QED result for the 
pions. However, once this point of view is accepted, the 
chiral enhancement looses its power as the theoretical 
parameter and the pion box contribution becomes just 
one of many 0{N^) contributions about which nothing 
is known at present. Therefore, for the final estimate of 
a'l^' we use 

a™'^° = 0(10) X 10-", (43) 
where the error estimate is clearly subjective. 

VI. CONCLUSIONS 

In this paper, we revisited the issue of the hadronic 
light-by-light scattering contribution to the muon anoma- 
lous magnetic moment, incorporating constraints from 
perturbative QCD in constructing the low-energy model 
for the light-by- light scattering. To achieve that, we com- 
puted the light-by-light scattering amplitude at a rela- 
tively large photon virtualities in perturbative QCD and 
required that low-energy models for hadronic light-by- 
light scattering have to interpolate smoothly between 
small and large values of q^. The minimal large- A'c model 
with such feature contains the pseudoscalar and the pseu- 
dovector meson exchanges. 

Since, by construction, the hadronic model we use in 
this paper has correct scaling at large values of g^, we 
achieve reasonable matching between the low-energy and 
the high-energy degrees of freedom that contribute to 
hadronic light-by-light scattering amplitude. One of the 
major findings in this paper is the fact that too strong 
a damping of hadronic amplitudes at large values of 
has been used in previous studies |E S 111 of hadronic 
light-by-light scattering to a^. 

It turns out that imposing correct matching between 
the low- and the high-energy degrees of freedom leads 
to substantial changes in both the pseudoscalar and the 
pseudovector contributions, making both of them larger. 
In a way, the impact of large- momentum degrees of free- 
dom on was underestimated in previous analyses. Our 
final result for hadronic light-by-light scattering contri- 
bution to is: 

a™ = 136(25) X 10"". (44) 



The error estimate includes the sum of 0{N^) error esti- 
mate in Ea. H43|l as well as 15 x 10"^^ as an error estimate 
for the sum of the pseudoscalar and the pseudovector ex- 
changes. From Eg . (|44|l it is clear that we do not claim 
significant reduction in the theoretical uncertainty, al- 
though we believe that it adequately reflects our current 
knowledge of aJJ^'. On the contrary, we think that the 
shift in the central value is real because it originates from 
a better matching of the low-energy hadronic models and 
the short-distance QCD. 

The result in Ea. l44|) is approximately 50 percent 
larger, than the currently accepted estimate ~ 80(30) x 
10"-^^ Note however, that our result is closer to 

another recent evaluation of hadronic light-by-light scat- 
tering contribution to a^^ [19| where the central value 
a™ = 108 X 10~" is quoted. 

Another possible consistency check is to estimate the 
light-by-light scattering contribution as a sum of two 
terms ~ the pion-pole contribution to account for low- 
momentum region and the massive quark box contribu- 
tion to account for large-momentum regime. If the quark 
masses are chosen to be ruq = 300 MeV, the result for 
the quark box contribution is 60 x 10"^^. Combining 
this with the pion pole contribution, we get the estimate 
ajjj'' « 120 X 10"^^. Of course, the above consideration 
is not the proof; yet it clearly indicates the tendency of 
the result for a!!'' to increase once the contribution of the 
large-momentum region is accounted for in the correct 
way. 

Finally, we note that the new value for hadronic light- 
by-light scattering contribution, Ea. (|l^ . brings the es- 
timate of the muon magnetic moment anomaly in the 
Standard Model and the current experimental value 
somewhat closer. Using recent results for hadronic vac- 
uum polarization we arrive at: 

oxp _ „th _ r (171 ± 110) X 10-" (1.5 a), 
M ""a" - I (24 ± 110) X 10-" (0.2 fx). 

The first result uses e+e--data only, while the second 
one uses the r data at low energies; the errors in each 
equation are combined in quadratures for compactness. 
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Appendix I 

In this Appendix we give explicit expressions for the 
light- by- light scattering amplitude in perturbative QCD 
in the kinematics when three photons have non-zero vir- 
tualities and one of the photons is soft. 

The Euclidean amplitude 

A — ,Ml M2 /i3 M4 A (Af^\ 
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defined in Eq.Q can be expressed through nineteen gauge- invariant structures: 
I 

A = G^''=^^{//i}{/2/3} + Gi'''''^{//2}{/3/l} + Gf''''^{//3}{/l/2} 

+ G^^^' [fh] [hh] + G^^'' [f~h] [hh] + Cr^'^ [fh] [h~h] 

+ G^'-'^ {7723/ /i m^} {f2h} + Gf {7731/ /2 mi} {hh} + cf'^^^ {^uf h V12} {hh} 

+ Gi'^'-'^ {m3f h qi}{h h} + G^3^'-'^ {mfh 92} {h h} + Gi'^'-^^ {7712/ h 93} {/i h} 

+ g'-^'^'^^ {7723/ h ms} {92/2 V3i} {qsh m^} + Gf {7/31/ /2 7731} {^3/3 7712} {qih ma} 

+ cf'^'"^^ {vi2f hvi2} {qihma} {q2hV3i} 

+ G^ '^'^^ {qif fi qi} {92/2 f73i} {qsh '712} + Gf {^2/ /2 <?2} {q-ih ^12} {qih V23} 

+ Cf'^'^^ {qsf fs qs} {qifi 7723} {(72/2 Vsi} + G^^'^^^^ {qif qs} {qifi 7723} {^3/3 '712} {^2/2 ?73l} • (47) 



We have introduced the field strength tensor for all of 
the four photons = gf ej- - ef with f^" = /^"^ 
and, also, the four- vectors ijij — qi ~ qj. In Ea. (|47|l . 
we view /[''^ as matrices; the curly brackets imply either 
traces of matrix products or convolutions with vectors 
qi,Vi]- For example, {91/92} = qi,f,f^''q2,^- The no- 
tations for invariant functions G1-5 are introduced for 
compactness. The superscripts denote the arguments of 

these functions, e.g. Gg^'^'"^^ = G5((7i , g|, The invari- 



ant function G5(gj^, qf, is totally symmetric with re- 
spect to the permutations of its arguments; the functions 
Gi,2,3,4(9i) 92i 93) ^I'S symmetric under the permutation 
of the last two arguments; the functions G3,4('Zi , f?! ' 93) 
are antisymmetric under the permutation of the last two 
arguments. 

We have computed the above form factors in pertur- 
bative QCD where the photon-photon interaction is me- 
diated by the loops of massless quarks. Our results are: 



Gi(si, S2, S3) = - (453^ - 4si^s^ + Gsjsl ~ ^szs\ - 22sjs^ + 2'^,sj,s\s\ + 2%s\s\sy ~ \\s\sx + 32s?s^ 
-48S3S2S1 - 2s2 -I- S3S1 -I- 2siS2 - 4S3S2S? - 32S3S2S1 -I- 32s|s2Si + 2s|s2 + Si - ^s\s 



3g2N ln(s3) 

^ ^' D'^Si{si - S3 - S2)S2 



+ {si - 4s^si + Gslsl + 6slsl + 2IS34 - 22.s^s^ - 22si^s^ - 26S3S2S? + 56slslsi - As^sl + 2IS3S2 

+4S3S2S? - 26sls2sl + si + sis2 - -islsl - 4^51 + Sgsf)— — —i^^ifi^ — 

-D S3(si - S3 - S2)S2 

- (-4sfs2 - 32s3s|s^ -I- 28S3S2S? + 32s3S^si -t- 26s^s^si - Asjs^ - 2s| -I- 6s?s| - 4s§sf -1- s|si -I- 2s3s| 
-h4s|s2 - IIS3S1 +sl- 4S3S2S? - 48sis2Si - 22s^s? -I- 32sj^s? + 2s3S^) 

13^5153(51 - S3 - 52) 

-2 (I8S3S2 -I- 75^51 - 7sjs2 + I8S3S251 - I953S2S1 - 45^52 "t" 35^53 - 4535^ - 75^53 + 2sj + isjsl 

3 2 4 . 4\ -^(S1,S2,53) (-5^ - 3siS3 5? - 5^ - 352Si -t- 2S352) 

+ 751S| - 19S35251 - 553 - 5S2) ^r27 T " 4 



D2(si-S3-52) (Sl - 53 - 52)siI3 

G2(5i, 52, S3) = -8 (252 ~ 453S2 + 253 - 5251 " S153 - 5^) ^-^^J^ + 4(s2 - Si - 53) (sg - 25352 
-25153 + 5^ -f 452S1 + sl) - 4(-S3 + 52 + Si) (sj - 2s^S2 + 4siS3 + 5^ + 5^ - 252Sl) 

U-^si U'^si 

0/ 3 2 2i3r,2ir, r,2i2i2 \ "^(Sl, S2, 53) 8 
-8(S2 - S253 - 5352 + 53 - 2si52 + 2515352 - 25i53 + S-^S2 + S1S3) — , 
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G3(si, S2, S3) = 2 (5.sJ|.si - lOslsj - 26slsl + 22s^S2S^ - 2OS3S2S? - sls2 + dislsl 
+10slsl + 15s^ - 474.S1 + BOsjsl - SSssslsf - ISsfsl + 20s^S2Si - 21s3S^ + 84s3sisi 

-62slslsi - 5s3st - 21^52 - 4 + si) 

-2 (Ss^si + S3sl - 6^52 - lOsjsl -4-4- 38s|.S2Si - SSsss^s? + 22s^s^si + 22s|s^si 
-27s3S^si - 2s3S2sf + 72sls2sl - 76s^s^s? + 72s3sf .s^ - 27.sj.S2Si + .sf S2 + SSs^.s^ 

-6S3S2 - 52s|s^ + 33S3S2 + 5slsi - IOS1S2 - ^44 - Ssgst + lOslsl + 10s?s^) ^^^^^ 

+2 (5s^si + 22s3s|sf - S2 - 62s^s^si + 20s3S^si + - 21s^S2 - 26s|s^ + 15s| + lOsfs^ + 84s;^.s2Si 

-38s^S2S? - lOsjsl + 54:slsl - 18s^sf + 50s^s? - 20s3S2sf - sss^ - 5s^S2 - 47s^si - 21s3S^) ^^^^ 

+4 (-8s?s^ + ISslsl - ASslsjsi - 4s3slsl - 9s3s| - 9^52 + 3s^ + 3s| - llsss^ - 4s^si + 2sf 

+28s3.s^si + 18s^s? - Sslsj + 28s|s2Si - 4s^si - Il4s2 + Gs^sl - 4s3S2S? - 4s^S2Si + 6s^s^) '^''^^^3'^'^'' 
8(sf - 4sl + 8S3S2 - 4si + 351S3 + 3s2Si) 

£)2 

G3(S1, S2, S3) = -2 (-S3S1 - 5S3sf - 2s|s2 + 5S3S1 - 2sl + S? - 3OS3S2S? + 26S3S2sf + 98S3S2S1 - IOS3S2S1 
-53S3S2S1 + IOS3S2S1 + 56S3S2S1 - II4S3S2S1 + I6S3S2S1 - 27S3S2S1 - 7sf S2 - 4S3S2 + 6S3S2 - 4S3S2 + 6S3S2 

-7s|si + 37sfs| + 28sts^ + 10s?sf - 10s?s? - 50s?s^) MfA^ _ 2(s2 - s^) (s^ + 5s^si - 2Ls3.s^ + 24535^51 
-lOsfsl + 22slsl - IO6S3S2S1 + lOsfsl + 10s3S2S^ + 22S3S2 - 2IS3S2 - 5s^S2 - 8s3S2S^^ + 24s^S2Si 
+ 10s^s2s2 + 5slsi -si- lOslsl + sl + lOsjsl - SsaSi )^^^^ + 2 (-7s|si - 7s3S^ + 6s|s2 + 37s^s? 
+sl - 2sl + 26s^S2S? - 3OS3.S2S? - lOslslsi + OS^SaSi - 27S3S2S1 + IOS3S2S1 + 16s|s2Si - II4S3S2S1 

+ 56.S3S2S1 - 53.S3.S2S1 - 5.S1.S2 + 6.s|.S2 - 2.S3.S^ - 4s|s2 - 4S3S2 - S2S1 

+5s?s| + lOs^s^ + 28slst - 50s|s? - lOsfs^) j^^^ - 4(s2 - S3) (-5s^ + 14s^si - 4s3S^ - 12s^s^ 
-26s3sisi + I8S3S2 + 28S3S2S1 + 2sfs2 - 4s^S2 - 26S3S2S1 + Si - I2S1S3 + 14sis^ - Ssg + 25^53) 

8(S2 - S3) 



(-S2 + 2S3S2 - 3S2Sl - S3 + Asj - 3siS3)si£)2 ' 



G4(S1, S2, S3) = -4 (-4 - 10S?S2 + 8S3S^ - 2s?S3 + 2sisl + Si - 17S2 + 48S3S2Si - 78S3S2S1 

I + 2s; 
In(si) 



-14s^S2Si - 24s|s2 + 26sisi + 34s^si) — ^-^^^ -— - 4 (-s| - Ssgs^ + 2sisi + 18s^si 

S2(si - .S3 - .S2)(.S2 - .S3)£''=' 



-38s3sisi - 8S3S2 + 48S3S2S1 - 2s?S2 - 38s?S2Si - Sq + Si - 2sf S3 + 2S1S3) ^ , , 

L'-^S3(Si - S3 - S2)S2 

+4 [si - 78S3S2S1 - S2 - 2s?S2 + 2s^si - 10sfs3 + 26sis^ - 17s3 + 34s^si - 24s3S^ + 8s^S2 

+48S3S2S1 - I4S3S2S1) 7 Mf2) ^ _^ 24 (si - S3sl + S2S1 + IOS3S2S1 - 5SiS2 - S3S2 

(Sl - S3 - S2)(S2 - S3)S3i^'^ 

,2 , 3 , q3 r 2\ J{si, 32,83) 8(si - 2s2Si + 8S3S2 + sf + S§ - 2S1S3) 

+S3S1 + S3 + >5Si — OS3S1 1 —r- r -\ -, r , 

' D'^{S1 - S3 - S2) £'2S3(S1 - S3 - S2)S2 



G^isi, S2, S3) = 4 (-2S3S1 - 22slsl + Asjsl + 4s2Si - Asssf + 22S3S2S1 - 26S3S2S? + 48s3S2Si - 2S3S2 
-26s3slsl + IIS1S2 + 2S3S2 + 4S3S1 + 2sf - 32S3S2S1 - S3S2 + S2 + 8S3S2S1) 



Sl(si - S2 + S3)D^ 



2 



-4 (2S3S1 + Asjsi + 26S3S2S1 + 2S2S1 - S2 - 4s3sf - 2s| - 48s^S2Si + 26S3S2S? - 22S3S2S1 - 8S3S2S1 

+22slsl - Il4s2 - 2slsl + 32s3S2sl - Asjsl ~ As^s^ + sfs2) ^^^7^3 2 

+4 (-2siS2 - 2S3S2 + 52S3S2S? + sl + 2IS3S1 + 2IS3S1 + 2S3S2 - S3S2 - S2S1 - 8S3S2S1 - 22slsl 
-22slsl + 2slsl - 12s3slsi - 12slslsi + s|) ^""^^"j ,^3 

S2Sl(Sl - S2 + SsjD'^Ss 

+8 (l7s|s2Si + I7S3S2S? - 5S3S2 - 5sf S2 + 2s| - 18sf Sg + S3S2 + Ssg + 4sf S3 + 4siS3 + 5sf 

n 2 2 q22|3 in 2\ J'(S1,S2,S3) 

+ 



S2(.'5l - S2 + S3)I?3 

S2S1 — 2s^S2 + ■s?S2 + S3S2 — 2s§S2 + S3S2 + 2S1S3 — 4s^s| + 2s5s3 + 5s|s2Si + 5s3S2sf — 2s3s|si) 



sls3D^{Si - S2 + S3)si 



G5(Si,S2,S3) =7?.i(si,S2,S3) + T?-! (S3, S2, Sl) + T^-l (si , S3 , S2) + 7?.2 (Si , S2 , S3) ; 

7^l(sl, S2, S3) = 8 (s? + s| + s^ - IGslsi - 513sf S3S^ - 181s^sf + 1737sf s^s^ 
-515sf s^S2 - 16s|s| - 10s^S2 - sgsf + 26s^s^ + 26s|s^ - ISsgS^ - s|s2 - sssf + 26s|sf - lOsgS^ 
-10s|sf + 26s^sf - 16sfsf - lOs^s^ - sis| - 125s?s|s2 - 1544s?s^s| - 453sfs|s2 - 446s?S3s| 
+1862s?s|s^ + 1737s?s^s^ - 990s?s^s| - 1645s?s^sf + 177sis^S3 - 183sis|s^ + 207siS2s| 
+909sis|s^ - 125sisisf - 515sis^si - 453sis^s| + 251s?s| - 55s?s^ - 16sis| - 181sts^ - SSs^s^ 
+251sfs^ + I566S1S3S2 - 990sis^si + 909S1S3S2 - 1544sts3si - 1645s?S3si - 513s?S3S^ + 3550s?s^s^ 

-16sfs2 + 177s[s3S2 - 446sfs3S^ - 183sfs^S2) ^^^^^"^ 



D^DiSlS2{si - S3)(S2 - S3) 

^2(si, S2, S3) = -16 (-54s|si - 54s3sf - 54s|s2 - 45s^s? + 29s| + 29sf + 29s| - 288s^S2S? 

-288s3S^sf - 288s^s^si - 288s^s^si - 288s|s^si + 342s3s|si + 342s3S2s| - 288s^S2sf 

+666S3S2S? - 288s3sisi + 342S3S2S1 - 54sf S2 - 45s3s| - 54s3S^ + 140s|s^ - 45s3si - 54s^si - 45sis| 

-45s^2 _ 45sisf + UOslsl + UOsfsl) 

-16 (sf + 58s3S2sf - 9sfs^ - 9s^s^ + 16s^sf - 58s3s|sf + 16sfs^ - 58s^S2sf - 9s^sf 
+134s3s|sf - 58s3S2sf - 9sf S2 - 58S3S2S1 + 58s3s|si + 58S3S2S1 - 58S3S2S1 - 58s3s|si + S3 

-9S3S2 - 9S3S2 + S2 + I6S3S2) £,3£)^3^32S3' 

In the formulas above 

D = S3 + Sl + S2 - 2siS2 - 2siS3 - 2S2S3 , 

-Dl = (S2 + Sl - S3)(S1 - S3 - S2)(S1 + S3 - S2), 

and the function J(si, S2, S3) is defined through 

d^Z 1 



J(si 



,S2,S3) = j 
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where = Si. The function J(si, 32,33) is symmetric w.r.t. to all its arguments. The explicit expression for this 
function in terms of the polylogarithms of rank two can be found in pol |. 



Appendix II 

Below we give the results for the functions fi{5, L) for i = 0, 1, 2 introduced in Ea. (|38|l : 



11 16 



11 



1925 



/o((5,L) = -— - —04 + — CsTT^ " -Cs - -Cs + 127r2ln2 + -TT^ In^ 2 - 



+5 



+5^ 



36 



72 



1604 - 4C3^' + 6C3 + |C5 - IStt^ In 2 - ^ 
18 6 3 



943 2 

108" "540" 



31 4 2 4„ 
TT In 2 

540 9 



- ' ?ln^2 



64 5 , 
3 12^"^ 



2 l^A 

24 ^' 



5. 47 „ 
3-tC5 + ^^' ln2 



9 



TT^ In^ 2 - 



54 



540 



ln4 2 



. ,7 161 104 2, „ 



_ i5_^4 24a4 - TT^ In^ 2 + In^ 2 - -Cstt^ 
648 216 9^ 



^C5 



(51) 



/i('5,i) 
5^2 



13 27r2\ 29 40 27 4C37r2 67C3 2OC5 

LL H 04 ^2 H 

4 3/ 93 8 3 6 3 



9 

37 



ln^2- ^^"^ 



3l7r'^ 



27 



C3 



+ 



216 9 



10 21 2^^ 127 2 

TT In 2H TT 

9 54 

7^2 



In-* 2 + <5 [ 3^2 



2 2 1\ , 97 80 27^ 
' — — —a^ -02 



11 



-TT 



l^ln^2 
9 



^^^^^ 7r2ln2 



24 9 



\n 2 



60 

537r2 IStt" 



36 
1 

-TT 



2 3\ ^ 115 56 27 



^2 



72 



216 



-In" 2 
9 



329^0 /2597r2 



36 



V 72 



6C3 



14813 27 \ 40915 
432 ~ y V ~ 1728 



16a4 —S2 —^2 



1547 2 2, 2^ 217 2 

Cs In 2 TT 

36 3 162 



^." + ^^n" / 110^3 _ 130^2 

4320 3 V 9 9 V 36 



1257r2 5239 



108 



27 \ 30175 

Jo I J-J — 

4 / 1296 



32a4 — ^S'2 g~'^2 



7r3 673, 4 2, 2^ 160 2 101 4 4, 

^ + C3 + --^ In 2 H TT^ H TT* - - In* 2 

18\/3 18 3 81 216 3 



(52) 



(53) 



Here, L — ln(Afp/rn^), S — {rrif^ — m^)/m7r, C™ ^'^^ the Riemann zeta functions, 04 = Li4(l/2) and 5*2 = 
0.260434137632161. I 
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